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The interaction of a gravitational wave with a system made of an RLC circuit forming one 
end of a mechanical harmonic oscillator is investigated. We show that, in some configurations, the 
coherent interaction of the wave with both the mechanical oscillator and the RLC circuit gives rise 
to a mechanical quality factor increase of the electromagnetic signal. When this system is used as 
an amplifier of gravitational periodic signals in the frequency range 50 — 1000 Hz, at ultracryogenic 
temperatures and for sufficiently long integration times (up to 4 months), a sensitivity of lO"'^*- 
10"'^^ on the amplitude of the metric could be achieved when thermal noise, shot noise and amplifier 
back-action are considered. 
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I. INTRODUCTION 

Gravitational wave detection is a research field which is actively under development by many groups all over the 
world (e.g., From the theoretical point of view all the devices yet operating or under construction could be 

divided into two classes: the resonant detectors and the interferometers. In the first case one tries to detect the 
motion of an elastic body due to the gravitational tidal forces acting on it. In the second one the gravitational 
phase shift on an electromagnetic wave is exploited. In both cases the detection is due to only one kind of direct 
gravitational interaction with matter. For instance, in the resonant bar detector arrangement, the direct interaction 
of the gravitational wave with the solid body produces mechanical oscillations which are converted into an electric 
signal by a capacitive transducer and then amplified and recorded. The bar is the only constituent of the apparatus 
which directly interacts with the gravitational wave. The transducer system can not distinguish if the mechanical 
motion is due to a gravitational wave or another source. 

In this paper we will study the behaviour of a mechanical oscillator and an electromagnetic circuit rigidly mounted 
on one end of it, both directly interacting with a periodic gravitational wave (see Fig. 1). As a first preliminary result, 
we will find the rest frame of one end of a damped oscillator in the field of a gravitational wave, generalizing a result 
contained in Then we will show that the device described in Fig. 1 is characterized by a sort of coherent effect: 
in some configuration the effect is to amplify the response of the electromagnetic system by the mechanical quality 
factor of the oscillator. This coherent effect, besides being interesting in itself, could be taken into consideration in 
future experiments for almost two reasons: in fact we will show that the detection limit of such a device, for periodic 
sources like pulsars (frequency range 50— 1000 Hz), at ultracryogenic temperatures and for sufficiently long integration 
times (up to 4 months) seems very interesting, being of the order of 10~^'^-10~^^ in the metric perturbation (when 
thermal noise, shot noise and amplifier back-action are taken into account); moreover, and this is the most important 
peculiarity, there is a decrease of the problem of mechanical insulation, because the sensitivity of the circuit can be 
considered independently of noise affecting the mechanical oscillator. This is due to the fact that in this arrangement 
the electromagnetic circuit, in some sense which will be clearer afterwards, "recognizes" if the mechanical oscillator 
signal is a gravitational one or not. 

II. GRAVITATIONAL WAVE INTERACTION WITH ELECTROMAGNETIC CIRCUITS 

Recently it has been shown ^ how, under certain conditions, a set of conductors behaves in the field of a gravita- 
tional wave. These conditions are essentially three. The first one is about the electromagnetic field: it is considered 
in the long-wavelength approximation in which the frequency of the electromagnetic oscillations z^em is much smaller 
than c/d {d is the typical linear size of the system and c is the speed of light). The second one is about the gravitational 
field: the spatial distance Xg over which the gravitational wave changes appreciably must be much greater than d. 
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The last condition is about the mechanical behaviour of the system: it is freely falling in the field of the incoming 
wave. This implies that it is at rest in TT-gauge, reference frame in which the calculation is performed. Within 
the framework of these assumptions the equations for a system of N conducting non-ferromagnetic one dimensional 
(wires) or extended bodies (capacitors) in the presence of a gravitational wave are derived. This is made by means of 
an extension of the lagrangian formalism (see, for instance, |^), in which the charges on the conducting bodies play 
the role of generalized coordinates. We have found that the behaviour of such a system is completely described by a 
lagrangian together with a dissipation function in which the coefhcients of mutual inductance L-y^^^, the inverse of the 
capacity coefficients C~\^^ and the resistances R^^ (t) that characterize the system are modified by the gravitational 
wave (here the subscript 7^ means the a-th conductor). This could be summarized by the following relations 

^lalb (0 ~ ^lalb ~^ '^lalb (^) 



(2.1) 



R-y^ {t) — ^R^ 



At) 



where the zero quantities are the usual flat space-time ones, while l-y^'y,,(t), c^^~i^{t) and T^^{t) are the time dependent 
perturbations induced by the gravitational waves on the circuit parameters; they can be written as 

^1alb{t) — ^ij (^)-^7„7i, 

c^.^At) = h^,ml^b (2-2) 

where hij{t) is the metric perturbation which describe the gravitational wave in TT-gauge (5^1/ is the metric tensor, 
greek indices run from to 3; latin indices from 1 to 3) 

r]f,^ = diag{-l, 1,1,1); \hf,^\ « 1 (2.3) 
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The coefficients AJ^^^^^, X^'-^^ib ^^"^ ^la 
role as the mass-quadrupole tensor 



defined in Eqs. (2.4)-(pT6|), are pure geometric quantities which play a similar 



p{x){x'x^ - {l/2,)5'^x''xk) d^x 



in the mechanical interaction between a gravitational wave and a solid body. We see therefore that the equations 
that describe the system are parametric ones. This means that gravity is not acting as an electromotive force but 
is changing the circuit parameters. The effect of the small perturbations of the circuit parameters when there is no 
electrostatic charge distribution is a frequency and amplitude modulation (e.g., |^). 

Let us now consider the circuit described in Fig. 2. It is an RLC circuit with two plane and statically charged 
condensers set at right angles. When a periodic gravitational wave (for the sake of simplicity coming from a direction 
perpendicular to the circuit plan e) in terac ts w ith such a circuit, a current is produced. In this case the effects of the 

and 



parametric interaction given by (2.1 



|2^) could be summarized by the relations 
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where L, C and R are the unperturbed inductance, capacitance and resistance and the perturbations are very small 
compared to unity, being of the order of the amplitude of the perturbation of the flat space-time metric in TT-gauge. 
In a reference frame in which the x and y axis are parallel to the plates of Ci and C2 respectively then the current 
flowing is a small modulation of the solution of the following equation (Ci = aC, C2 ~ C, = [{a + l)/a](l/LC), 
Qi = aQ2, Q2 = Qo/(a+l)) i 

Q + 2-iQ + Lol^^Q = v+h+[t) (2.8) 

where 

= (^^'^oOo (2.9) 



and hj^{t) is the "plus" TT-gauge polarization amplitude of a gravitational wave ||6[. Eq. (2^) is valid both for the 
ohmic and the superconducting case, provided that 27 is the inverse of the damping time of the circuit in the two 
cases. The conclusion is that in the case of interaction with circuits having a distribution of electrostatic charge the 
most important effect of the gravitational wave is to produce an electromotive force having the order of magnitude 
of the electrostatic potential in one condenser times the amplitude of the metric perturbation in TT-gauge. If we 



consider Eq. (2.9) we can see that it is the analogue of the tidal acceleration in the equation of motion of a mechanical 
oscillator. Here the electrostatic charge has taken the place of the rest length of the oscillator, and the capacitance 
that of the inverse of the elastic constant of the spring. 

In which follows we will study what happens if we rigidly mount with no moving parts this particular circuit on one 
end of a mechanical harmonic oscillator. The key idea is to treat the electromagnetic system in the frame in which it 
is at rest so we can use the results of Ref. H . 



III. THE REST FRAME 

In the local inertial frame x^, the so-called Fermi Normal Coordinates (FNC frame) 1^,0, the interaction with a 
gravitational wave of a harmonic oscillator made of two masses M , connected with a spring of equilibrium length 2L 
in the direction {U = Ln'), in the case in which the wave depends only on time (L << \g), is described by the 
geodesic deviation law (see, for instance, 1^, p. 1023) 

f + = (3.1) 

To 2 

where (L* + C*) are the FNC of one of the two masses, wo is the proper frequency of the oscillator and tq its relaxation 
time. The right-hand side of this equation contains the "electric" components of the Riemann tensor of the wave in 
TT coordinates, R^a^j = {l/2c^)hij. In linear approximation the Riemann tensor is gauge invariant, i.e. it does not 

depend on the chosen coordinate system. The solution in Fourier space [^(w) = (27r)^^/^ J^°° ^(t) exp(— iwt) dt], in 
linear approximation, is given by 
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C{^) = --L=h^,{^) ^ , (3-2) 



Within the framework of long wavelength approximation [L <^ Xg), the FNC metric could be written as (see, for 
instance, p. 332) 77^,^ -I- 0{x^x^ /Xg) where |a;*| < L. Therefore up to {L/XgY terms the FNC metric could be 
safely taken as rj^^i, in the region of our interest. 

Now we introduce a new coordinate system which follows the motion of the mass. To this aim, we consider the 
following coordinate transformations 

C+00 2 

a;* = - (87r)-i/2yJ / du eyi^{iojt)^ i-h\{uj) (3.3) 

J -00 u^-u^ + iu/tq J 

xa = yo + ^(87^)-l/2(y^yJ/2c) / duj exp{iu;t)— ———h,j{uj) (3.4) 

The new metric (the rest metric) can be written, neglecting [d/XgY terms, as [see, for instance, Eq. (83.6) or 
Eq. (4.3)] 
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Vfj^^ + = rj^^ = / dwexp(icji)^ :-———h^„{uj) (3.5) 



In this frame, in the same approximation, the mass under consideration is at rest; in fact, caUing (L* + i]^) its 
coordinates, then geodesic equation becomes 

^Irf = (3.6) 

To 



We note that, within the range of the approximations made, the rest frame has the same properties (2.3) of a TT- 
gauge one. We furt her observe that if we consider the hmit tjQ and tq — > oo (free masses) then H^j^y = h^y and 
Eqs. ( |3.3|) and (3.4) become the transformation laws from the TT-gauge to the FNC to the first order in the L/\g 
ratio Hnl; in other words the rest frame becomes the TT-gauge (as it must be because it is the rest frame for free 
bodies). 

Let us consider now the case of an electromagnetic circuit rigidly mounted on one end of a mechanical harmonic 
oscillator. In the rest frame we are within the framework of the theory summarized in section 2. We can therefore 
apply the results of 1^ by simply substituting the /ly metric with the rest metric. Therefore the rest frame solves the 
problem of the simultaneous interaction of a gravitational wave with a mechanical oscillator and an electromagnetic 
circuit, being the frame in which only the coupling with the latter occurs. 

It is to be remarked that we are dealing with currents in electromagnetic circuits in TT-gauge. In actual experiments 
currents are measured in the laboratory reference frame (FNC). However, because of the charge conservation we can 
write 

where /^^C" jg |.jjg current measured in the FNC and I^^ the one measured in a TT-gauge (in this case the rest 
frame). Because of eq. (p.4[), the relation between currents in the two gauges reads 



tFNC 



I^^ 7-00 UjI-UJ^+ iuj/TQ 

cdx^J_^ ' ^ 'uj^-u^ + iuj/to ^ ' 

In the long wavelength approximation both first order corrections are negligible so that 

jFNC _ jTT 

(see also (H). 

It is important to observe that the situation envisaged in this section is not equivalent to the case of a mechanical 
oscillator driven by an artificial force. In fact the equation of motion of the oscillator in this case is given by 

f + ^+L.o¥ = r(0 (3.7) 
To 

where f^{t) are three functions of time. If we now introduce a rest frame, that is to say a frame comoving with the 
mass element of the mechanical oscillator over which is mounted the electromagnetic device, and calcul ate the metric 
tensor in this frame, we find thatthe metric is unchanged. In fact, if ^* = C{t) is the solution of Eq. (3.7), then the 
coordinate system ij*^ of the rest frame is given by 

x' = y' + C{t); x°^f (3.8) 

The new metric will then be given by (it has to be remembered that in the long wavelength approximation the FNC 
perturbation of the metric is negligible) 

h^j = 0; /loi = Ci,o; ^00 = (3.9) 
In the same approximation as before Ci,o could be neglected; therefore we are led to the conclusion that 

V = (3-10) 

This result shows the difference between an external driving force on the mechanical oscillator and a gravitational 
wave, difference which lies in the tidal nature of the latter. 
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IV. MONOCROMATIC WAVES AND COHERENT INTERACTION 



Let us consider now the case of a monochromatic gravitational wave of angular frequency ujg. The coordinate 
system y'^ can be now written as [r — (LUg/ajo): Q„i = {loqTq)/2 is the mechanical quality factor]: 



while the rest metric takes the form 



Hij - 2 



Bir, Q„) 



where 



A{r, Q,„) = - 



1 r^(l-r^) 
2(l-r2)2+r2/4Q^ 



4Q™ (l-r2)2+r2/4Q2 



(4.1) 



(4.2) 

(4.3) 
(4.4) 



In the resonant case, ujg — loq, Ail, Qm) — while i?(l, Qm) — Qm so that the rest metric at mechanical resonance 
becomes: 



rjimr) _ 2Qm ; 



(4.5) 



Therefore the interaction of the gravitational wave at mechanical resonance produces, on the electromagnetic circuit, 
a sort of effective metric perturbation larger by a factor 2Qm and 7r/2 phase shifted. 



V. INDEPENDENCE OF CIRCUIT SENSITIVITY ON NOISE AFFECTING THE MECHANICAL 

OSCILLATOR 

The amplification mechanism described in previous Section concerns only gravitational signals that coherently 
interact with both the mechanical oscillator and the electromagnetic device. In fact, if the mechanical oscillator is 
accelerated, the order of magnitude of the inertial electromotive force induced on the electromagnetic system is 

£in ~ —ad (5.1) 

e 

where is the electronic mass e its charge and a the magnitude of the acceleration and d the linear size of the 
circuit.. If the energy of the motion (for instance thermal noise) is comparable with the one due to the gravitational 
wave interaction then a — aj^Lh where h is the order of magnitude of the gravitational wave amplitude. 

Let us consider now the circuit of Fig. 2. Th e order of magnitude of the electromotive force induced by a 
gravitational wave with h amplitude is [see Eq. ( ^.9[ ) or [^] 

V « Voh (5.2) 

where Vq is the order of magnitude of the electrostatic potential in the condensers. The ratio between these two 
electromotive forces is 

2 2 

!!!£^^iLd^2x 10-1°^ (5.3) 
V e Vo Vq 

where we have put Ld — 1 m2, i^g is expressed in Hz and Vq in Volts. Even setting Vq = 1 V and i^g — 10'^ Hz we obtain 
— 10^^ for the ratio. In conclusion, in order to calculate the sensitivity of the device it suffices to consider noise only 
on the electromagnetic circuit. Hence, at mechanical resonance, the sensitivity in h is increased by the mechanical 



quality factor (see Eq. (4.5)] 



5 



VI. DISCUSSION 



In this Section noise is evaluated in order to obtain the sensitivity of such a device. Taking into account discussion 
made in Sec. V, only noise on the electromagnetic circuit is considered. We take into account thermal mechanical 
noise, thermal electric noise, back action of the amplifier and shot noise of the difference of the leakage current of the 
condensers. 

Mechanical thermal noise describes the relative variation of length of inductance (or distance between plates of 
condensers) caused by thermal noise. In the free falling mechanical approximation (see 0]), the amplitude of th e 



gravitational wave having a signal to mechanical noise ratio equal to 1 after an observation time tabs, is [see Eq. (4.5)] 




where fcs is the Boltzman constant, T the temperature of mechanical degree of freedom, and M , I and r are the mass, 
length and damping time of mechanical oscillators describing condensers or inductance from the mechanical point of 
view. The analogue for the electric thermal noise at electric resonance {ujg — ujem) is [see Eq. (6.44)]: 



hetn^^— / (6 2) 

QmVo y QemCojgtobs 

where Qem is the electric quality factor of the circuit, and C = Ci = C2. 

Another noise source to take into account for is the shot noise of the difference, AJ, of the leakage current of Ci and 
C2. The gravitational signal with a signal to leakage noise ratio equal to 1 after an observation time t^bs is [see pO| ] 
and §, Eqs. (6.15) and (6.26)] 



1 2eAI 

/i/n = 77^ Tr — \ —: (6-3) 

4 Sdm Vo G <dem V tobs 

where e is the absolute value of the electronic charge and we can set AI — 10 p A [lO| . 

Once the bar and the electromagnetic circuit are isolated from the outside world by means of mechanical filters and 
Faraday cages a further noise source is the back-action of the amplifier. However, at cryogenic temperature and for 
the best amplifiers, this is at most of the same order of magnitude of thermal noise. Therefore, the amplitude of the 
gravitational radiation having a signal to noise ratio equal to 1 is given by 



2 (hintn + Ktn) + (6.4) 



In order to obtain the order of magnitude of the parameters which enter in Eqs. ( |6.lD -( |6.4| ) had to be set. 
Attention must be paid to the fact that weight and size of the electric circuit are constrained by the mechanical 
harmonic oscillator to one end of which the circuit is rigidly mounted. A possible choice is: Qm = 10^, M = IKg, 
T = Isec, I = 10~^m, tobs = 10^ sec ~ 4 months, Vb — 10 kV, C — 1 /iF, Qem — 10'^ in the ohmic case and 10^ for 
superconducting circuit. Moreover we assume to operate at cryogenic (4K) and ultracryogenic (40 mK) temperatures 
and for gravitational wave frequencies of 60 and 600 Hz. The reason of this last choice is twofold: first, the possibility 
to realize resonant mechanical detectors operating at these frequencies ^Jljl; second, the presence of possible almost 
monochromatic sources such as crab pulsar [^ or millisecond pulsars [j^. In this range it is found that the leakage 
noise is negligibly small. Moreover for superconducting circuit, mechanical thermal noise is always the biggest noise 
source, while in the ohmic case, at the higher gravitational frequency, electric thermal noise becomes more relevant 
than mechanical one. The results can be found in Tab. ^. 

In order to make comparisons with the predicted sensitivities of ground based interferometric detectors such as 
VIRGO and LIGO (e.g., [|l|), for gravitational signals from pulsars, we introduce the so called characteristic amplitude 
he defined in Ref. [Q, Eq. (50); setting Aj^ the "plus" TT-gauge amplitude polarization we have: 

where i is the angle between the rotation axis and the line of sight to the source. Therefore the amplitude he of the 
weakest source detectable with a signal to noise ratio equal to 3 in 4 months = 10^ sec integration time, with known 
frequency, is therefore given by 
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V 5 1 + cos^ I 



The results for ohmic and superconducting circuit are given in Tabs. |I^ and II] . The incertitude takes into account 
for different possible orientations of the rotation axis. 

A comparison with LIGO sensitivities (see, for instance, |15| ) shows that, near 60 Hz, the expected sensitivity is of 
the order of he ~ 10^^^ for initial detector and he ~ 2 x 10^^ for advanced detector, while near 600 Hz the expected 
sensitivities are 3 x 10~^^ and 10^^^ for initial and advanced detector respectively. 

Summarizing, at ultracryogenic temperatures, a device such as that described in this paper could have the same 
sensitivity as a ground based interferometer of the first generation near 60 Hz, while for higher frequencies the sensi- 
tivity seems to be comparable also with that of an advanced interferometer. These preliminary results are encouraging 
and show that a device exploiting coherent interaction amplification could be interesting in the field of gravitational 
wave detection. 



VII. CONCLUSIONS 



We have found the rest frame for one end of a damped harmonic oscillator in the field of a gravitational wave. 
This rest frame has been used to obtain the current flowing in an RLC circuit with an electrostatic charge, rigidly 
mounted on that end and interacting itself with the wave. The result on the circuit behaviour could be described as a 
coherent effect between two electromotive forces always in phase: the inertial one (produced by the elastic motion of 
the harmonic oscillator in the field of the gravitational wave) and the one induced directly on the circuit by the wave. 
Together with this it must be noted that, not having a definite phase, inertial electromotive forces due to noise bar 
vibrations do not couple coherently with the circuit. Therefore an RLC circuit mounted on one end of a mechanical 
bar detector operating on the same frequency could be an interesting amplifier of gravitational periodic signals, as 
for instance those produced by pulsars and in particular millisecond pulsars |13|. This new way of amplify signals 
could reduce the problems of the mechanical insulation of the bar, because noise bar vibrations do not couple with 
the circuit with the same strength as the case of vibration of gravitational origin. 
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FIG. 1. An electromagnetic system is rigidly mounted on a side of a mechanical harmonic oscillator with proper frequency 
uiQ and relaxation time tq. Both of them interact with a gravitational wave. 



FIG. 2. RLC circuit. A static charge (Qi ■ Q2 > 0) is distributed on the plates of the condensers Ci and C2. 



TABLE I. Gravitational amplitudes with signal to noise ratio equal to 1 in 10 sec ~ 4 months integration time for different 
operating temperatures and gravitational frequencies, hmtn, h^tn and refer to mechanical thermal noise, electric thermal 
noise and leakage noise, respectively. 



Ohmic circuit 


temperature (K) 


gravitational frequency (Hz) 


hmtn 


hetn 


hln 


4 


60 


3.3 ■ lO"""* 


3.8 • 10""'' 


3.8 ■ lO""'' 


4 


600 


3.3 ■ 10""'' 


1.2 • 10""" 


3.8 • 10"^" 


0.04 


60 


3.3 ■ 10""'' 


3.8 • 10""" 


3.8 • 10""'-' 


0.04 


600 


3.3 ■ 10""' 


1.2 ■ 10""'" 


3.8 ■ 10"''" 


Superconducting circuit 


temperature (K) 


gravitational frequency (Hz) 


hmtn 


hetn 


hln 


4 


60 


3.3 ■ 10""" 


1.2 • 10""" 


3.8 ■ lO"''" 


4 


600 


3.3 ■ 10""" 


3.8 • 10""' 


3.8 • 10"^'' 


0.04 


60 


3.3 ■ 10""=^ 


1.2 • 10""' 


3.8 ■ lO"""" 


0.04 


600 


3.3 ■ 10""' 


3.8 • 10""** 


3.8 ■ lO"'''' 



TABLE II. Characteristic amplitude he of the weakest monochromatic signal from a pulsar detectable with a signal to noise 
ratio equal to 3 in 4 months integration time with ohmic circuit. The incertitude takes into account for different possible 
orientations of the rotation axis of the pulsar. 





lyg = 60 Hz 


Vg = 600 Hz 


4 K 


(2.6-5.2) • 10""* 


(1.0 - 1.9) • 10""" 


40 mK 


(2.6 - 5.2) ■ 10""" 


(1.0 - 1.9) ■ 10""" 


TABLE III. The same as Tab. but for superconducting circuit. 




lyg = 60 Hz 


lyg = 600 Hz 


4 K 


(2.6-5.1) • 10"""* 


(2.6 - 5.2) • 10""" 


40 mK 


(2.6-5.1) ■ 10""" 


(2.6-5.2) ■ 10""' 
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Fig. 1 
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Fig. 2 



